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ABSTRACT 

We present 7-180/im photometry of a sample of hyperluminous infrared galax- 
ies (HyLIGs) obtained with the photometer and camera mounted on the Infrared 
Space Observatory. We have used radiative transfer models of obscured starbursts and 
dusty torii to model their spectral energy distributions (SEDs). We find that IRAS 
F00235+1024, IRAS F14218+3845 and IRAS F15307+3252 require a combination of 
starburst and AGN components to explain their mid to far-infrared emission, while for 
TXS0052+471 a dust torus AGN model alone is sufficient. For IRAS F00235+1024 and 
IRAS F14218-I-3845 the starburst component is the predominant contributor whereas 
for IRAS F15307-I-3252 the dust torus component dominates. The implied star forma- 
tion rates (SFR) for these three sources estimated from their infrared luminosities are 
Mf,^aii > 3OOOM0?/r~^/i^g^ and are amongst the highest SFRs estimated to date. We 
also demonstrate that the well-known radio-FIR correlation extends into both higher 
radio and infrared power than previously investigated. The relation for HyLIGs has a 
mean q value of 1.94. 

The results of this study imply that better sampling of the IR spectral energy 
distributions of HyLIGs may reveal that both AGN and starburst components are 
required to explain all the emission from the NIR to the sub-millimetre. 
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1 INTRODUCTION 

Since the discovery of the 'hyperluminou s 



IRAS P09104-f 4109 (iKleinmann et al. 19881) and IRAS 
F10214+4724 ( [Rowan-Robinson et al. 1991D the nature 
of these extremely luminous sources has been a matter 
of much interest. Objects with luminosities in excess of 



* Based upon observations with ISO, an ESA project with instru- 
ments funded by ESA Member States (especially the PI countries; 
Prance, Germany, the Netherlands and the United Kingdom) with 
the participation of ISAS and NASA 
f E-mail: verma@mpe.mpg.de, present affiliation MPE 



Lfir > 10^^'^^ h^Q Lq^ exceed the luminosities of normal 
spiral galaxies by three to four orders of magnitude and 
emit more than 95% of their bolometric power in the in- 
frared wavelength range. 

As with the lower luminosity ultraluminous infrared 
galaxies (hereafter ULIGs, Lfir > 10^^ L©), the underlying 
power source responsible for such tremendous luminosities 
is a topic of much debate and remains controversial (the 
so-called 'Starburst- AGN Controversy'). The space density 
of HyLIGs was found to be similar to that of quasars of 



^ as per th e definition for Hvpcr luminous Infrared Galaxies 
(HyLIGs) in |Rowan-Robinson (200o|). 
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compar able luminosity. This result led to the postulate that 



such sources are powered by deeply dust enshrouded AGN 
and may rep resent an evolutionar y stage in the formation of 
quasars (e.g. Banders et al. 1988a W) . Alternatively, the high 
infrared luminosity has been attributed to highly obscured 



compact nuclear starburst events (e.g. Soifer et al. 1984 



Rowan-Robinson et al. 1984 



Condon et al. 19911 ). IR AS d& 
tection s of optical and/or radio selected quasars (e.g. Miley 



Edolson & Malkan 1981: ; Rowan- Robinson & 



d 1989) revealed their IR emission to be dominated 



et al. ] 
Crawfo: 

by a broad component spanning 3-30/im. In contrast, the IR 
emission of starbursts w as found to be strongly peaked in the 



^60-10 0/jm range (e.g. Soifer et al. 1984; Rowan- Robinson 
et al. ] 984 ) . The spectral energy distributions (SEDs) of 



many luminous infrared galaxies (hereafter LIGs) display 
both peaks. The commonly accepted interpretation is that 
the IR luminosity is powered by a combination of the two 
mecha nisms. The presence of a deeplv dust enshrouded AGN 



(2OOOI , hereafter RR2000) estimate that there are 100-200 
HyLIGs over the whole sky with Seo > 200mJ2/|^. To date, 
only > 40 have been discovered. A number of HyLIGs have 
been found via correlations of the IRAS catalogues with op- 
tical, radio or known active galaxy catalogues leading to a 
significant bias in the sample of known IfyLIGs. Thus the 
high proportion of AGN-like sources in any HyLIG sample 
may simply be a selection effect. RR2000 compiled the first, 
large sample of HyLIGs selected from unbiased IR surveys. 
In a sample of twelve HyLIGs discovered from direct follow- 
up of IRAS or 850/im surveys, 50% have AGN-like optical 
spectra. In five of these cases, the MIR bolometric power 
originating from the AGN exceeds that of the starburst. The 
proportion of HyLIGs containing AGN is in agreement with 
the fraction (~ 5 0%) of Seyferts foun d in ULIG samples with 



Lfir > lO^'^-'' (Veilleux et al. 199S). However the fraction 



cannot )pe ruled out in any ULIG/HyLIG. However the pres- 
ence of an AGN does not imply that it powers all of the 
infrared emission. 

The role of ULIGs and HyLIGs in galaxy evolution 
remains unclear, as does the evolutionary connection be- 
tween the coexistent starburst and AGN activity detected 
in these luminous sources. ULIGs and HyLIGs have been 
proposed to be the dust enshrou ded stage in the evolution 
of a quasar dSanders et al. 1988a] Jb|). In this scenario a form- 
ing AGN is shielded from detection by an optically thick 
dust sc reen. As the AGN becomes more powerful the screen 



of sources with bolometrically dominant AGN i s higher than 
that found in ULIG samples (e.g. ~ 20 — 30% Genzel et al 



199q ). 

In this paper we use infrared photometry to investi- 
gate the nature of HyLIGs by comparing the resultant in- 
frared SEDs to theoretical models. In order to do so, the 
broad-band IR emission needs to be well measured at dis- 
crete wavelengths over the infrared range. However, for 
many HyLIGs emission in the MIR-FIR regime is poorly 
sampled from IRAS data alone, most sources only having 
one or two reliable]^ detections in the IRAS catalogues. We 
have therefore undertaken a photometry pr ogram spanning 



6.75-180/im using the camera {ISO-CAM, Cesarsky et al 



begins ]o quicklv break down revealing the optical quasar 1990 and imaging photometer {ISO-PHOT, 



{ Tanigi chi, Ikeuchi, fc Shioya 1999|). Alternatively, HyLIGs 
have b een interpreted within galaxy unification schemes as 



Lemke et al 



1996) instruments on board the Infrared Space Obs erva 



tory {ISO). We compare the data to pure star burst (Efs- 



'misdir ected' QSOs, the high luminosity analogues of Seyfert tathiou, Rowan- Robinson fc Siebenmorgen 200C| , hereafter 



2 galaxies (e.g. Barvainis et al. 1995; Hines et al. 1999| ). As 
yet, th e ubiquitv of AGNs in ULIGs/HvLIGs has not been 



conclus vely demonstrated although an increase in the num- 
ber of ULIGs displaying Sey fcrt-like characteristics has bee n 
recorded beyond L > 10^^'^ ( |Veilleux, Kim fc Sanders 1999| ). 

The resolution of the starburst-AGN controversy has 
implications for the star formation history (SFH) of the 
LIG population. In order to make accurate estimates of ob- 
scured star formation rates (SFR), the luminosity due to 
AGN and starburst activity needs to be quantitatively dif- 
ferentiated. The importance of this issue becomes apparent 
when one considers the popular belief that ULIGs/HyLIGs 
are the local analogues of the high redshift sources be- 



ERSOO), quasar ( Rowan- Robinson 1995, hereafter RR95) 
and inclined dust torus AGN~( pfetathiou, Hough, fc Young 



1995, hereafter EHY95) models to ascertain the nature of 



ing detected in sub-mi llimetre surveys (e.g. Hughes et al^ 
1996 ; f dson et al. 199^ ) . These dusty, high redshift sources 



may sample the populations contributing to the FIR-sub 
millimetre background. Such sources have also been linked 
to the X-ray background which, at hard X-ray e nergies, is 
postula ted t o have an origin in obscured AGN (Comastri 
et al. 1995; ^abian et al. 1998 ). Ultimately the resolution 
of the 'Starburst-AGN Controversy' may indicate a con- 
nection between these extragalactic backgrounds. Therefore 
ULIGs/HyLIGs may be exploited as local laboratories to 
investigate and understand the physical processes occurring 
at higher redshifts. 

Whether HyLIGs simply represent the high luminosity 
fraction of the well investigated population of ULIGs re- 
mains to be conclusively clarified. Yet, th is has proved dif- 
ficult due to the rarity of known HyLIGs. Rowan- Robinson 



the sources. Strong and broadly fiat MIR emission will in- 
dicate the presence of an obscured AGN. The contribution 
of the AGN component to the bolometric power can be im- 
mediately compared to that in the starburst component by 
using SEDs plotted in uFi,. 

We present photometric data for four objects from a 
sample of HyLIGs all of which have IRAS detections at 
60^m. These new data (a) enable detailed modelling of the 
broad-band emission for each HyLIG in this sample (b) pro- 
vide constraints for and allow the development of generic 
models for this class of infrared source (c) provide estimates 
of dust masses and star formation rates for these sources. 

The structure of this paper is as follows. Section 2 de- 
scribes the sample selection and is followed by an explana- 
tion of the data reduction methods in Section 3. Section 4 de- 
scribes background and models used in fitting these sources. 
The /SO- Cj4M images and spectral energy distributions are 
presented in Section 5. The results are discussed in Section 



van der Werf ot al. (1999) estimate a density of 7 X 10 deg 



for HyLIGs with Seo > 200mjy. This figure is higher than the 
estimate of Rowan-Robinson but is highly uncertain since the 
density is based upon the detection of only one HyLIG in an area 
of 1079 square degrees. 

i.e. flagged as good or moderate quaUty detections in the IRAS 
catalogues 
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6. We also address as sub-sections of the discussion (a) the 
implied star formation rates and (b) the extension (in terms 
of power) of the known FIR/radio luminosity correlation to 
HyLIGs by combining the IR luminosit ies with radio data 



taken from the NRAO VLA Sky Survey ( Condon et al. 1998 



hereafter NVSS). Finally, a summary of our conclusions is 
presented in Section 7. 

Unless stated otherwise, throughout this paper we 
adopt Ho = 50/i5o fems ~ ^ Mpc~ ^, Jlo = 1. 



2 ISO SAMPLE SELECTION 

The original ISO sample of HyLIGs contained all known 
HyLIGs (fifteen at that time) which were selected using a 
variety of methods and sampled a full range of optical and 
radio properties. Each object had at least a 60/im detection 
in the IRAS catalogues and the sample has a median red- 
shift of ~1. Due to scheduling constraints only five of these 
objects were observed^ The programs from which these five 
sources were discovered are described in the ensuing para- 
graphs. 



F14218+3845 is one of five HyLIGs found within this sur- 
vey which together efi'ectively represent a FIR flux limited 
sample of hyperluminous objects. 



2.4 IRAS F15307+3252 



IRAS F15307+3252 was discovered by putri et al. (1994| ) 
within a program investigating object s from the IRAS FSG 
'Warm Extragalactic Object' survey (Low et al. 1985) with 
faint optical counterparts. The spectrum of this source ex- 



hibits Seyfe rt 2 characteristics ( Kleinmann et al. 1988 
et al. 1994|). 



Cutri 



It is important to note that IRAS F15307+3252 and 
TXS0052+471 were discovered using selection techniques 
biased towards AGN (i.e. samples of galaxies with 'warm' 
IRAS colours and the cross-correlation of a radio- loud 
galaxy with the IRAS catalogues respectively). Whereas, 
IRAS F00235-H024, IRAS F14218-H3845 and IRAS F23569- 
0341 were discovered in systematic surveys, free from any 
bias towards AGN. 



2.1 IRAS F00235+1024, IRAS F23569-0341 

IRAS F00235-f 1024 and IRAS F23569-0341 were discovered 
in a systematic, optical identification program carried out 
with the APM machine at the Institute of Astronomy, Cam- 
bridge by McMahon et al. This project uses robust statistical 
estimators to identify IRAS sources taking into account the 
IRAS error ellipses and optical magnitudes of all potential 
optical counterparts. Deep VLA observations were taken of 
the sources with faint optical associations to confirm the cor- 
rect optical counterpart to be spectroscopically followed-up. 
These multi- wavelength investigations led to the discovery of 
IRAS F00235+1024 and IRAS F23569-0341, both of which 
display narrow line starburst spectra. 



2.2 TXS0052+471 

A cross-correlation of the entire IRAS Faint Source 



Data base with the TEXAS radio survey (Douglas et al 



1996) performed by Dey & van Breugel (1995), yielded a 



sample of radio loud IRAS galaxies of which TXS0052+471 



is one of t he most luminous [at a redshift of 1.935 (Jarvi 



et al. 2001 )]. Our ISO observations confirm the IRAS reject 
catalogue detection (FSR0052-I-4710). 



2.3 IRAS F14218+3845 

IRAS F14218-I-3845 was identified within the same ob- 
servational pro gram which led to the disco v ery of IRAS 
F1021 4+4724 ( aowan-Robinson et al. 199l| ). [Oliver et af 



(1996) carried out a system atic redshift survey of 1438 
IRAS Faint Source Survey ( iMoshir et al. 1992| ) sources 
with Seo >0.2Jy within a 720 square degree area. IRAS 



* A further five HyLIGs were observed under the 750 central 
program but are not discussed in this paper. 



3 ISO OBSERVATIONS AND DATA 
REDUCTION 

Observations were carried out using ISO's imaging photome- 
ter array ISO-PHOT at FIR wavelengths and imaged in the 
MIR with ISO-CAM. The majority of HyLIGs have two or 
fewer good or moderate quality IRAS detections. Thus ISO 
observations have been the only means, to date, of obtain- 
ing better sampling of the broad band IR emission of such 
sources. 

Images of the five galaxies were taken using ISO-CAM 
in CAMOl mode which was normally used for photomet- 
ric imaging. Observations were carried out using the LW2 
(centre 6.75/im, range 5.0-8.5/im) and LW3 (centre 15/im, 
range 12.0-18.0/im) filters with an integration time of ~300s. 
These images were used to (a) determine whether the ob- 
ject is extended in the NIR/MIR (b) obtain more accurate 
positions than IRAS (c) obtain accurate MIR fiuxes derived 
from aperture photometry from the images. 

ISO-PHOT observations were taken using the obser- 
vation modes PHT03 at 25/im and PHT22 at 60, 90 and 
180/im. These filters were chosen to span the greatest wave- 
length range and complement the existing IRAS data. 



3.1 ISO-CAM Imsiges 

Data reduction was performed from the raw data stage us- 
ing the CAM Interactive Analysis software^ (CIA). Firstly, 
the actual observations were extracted from the raw data 
product producing a data cube [image(x,y) against expo- 
sure time]. The data cube was processed with corrections 

^ The ISO-CAM data presented in this paper was analysed us- 
ing 'CIA', a joint development by the ESA Astrophysics Division 
and the /SO-C/IM Consortium. The /50-C^M consortium is led 
by the ISO-CAM PI, C. Cesarsky, Direction des Sciences de la 
Matiere, C.E.A., France. 
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applied for dark current, automated and interactive removal 
of cosmic rays and their remnants (deglitching) and sta- 
bilisation of the detector. The images were normalised to 
ADUs/gain/second. The data cube was then reduced to a 
single e xposure /image. This single image was flatfielded us- 



ties associated with the background model, the 25/^m flux 
was greatly over-estimated in comparison to the IRAS lim- 
its. As a check, IRAS and COBE/DIRBE background es- 
timates from IRSKY were also subtracted but ga ve simi- 
larly low backgrounds to those calculated from the Rowan- 



ing the internal CIA ISO-CAM flatfield. Finally, the CIA Robinson et al. (1991) model. Therefore, we concluded that 



structures were written out to a FITS image file with astro- 
metric calibration. 

The resolution of the KO- CAM images (6" x 6") was 
insufficient to resolve any structure or morphological infor- 
mation. Since the images were astrometrically calibrated, it 
was possible to determine the position of each source with 
an accuracy superior to that of IRAS. Aperture photometry 
was performed to obtain the NIR/MIR fiuxes using a fixed 
circular aperture with a diameter of 4 pixels^ centred upon 
the source. Surrounding this aperture, the background flux 
was determined using an annulus with an inner and outer 
radius of 2 and 6 pixels respectively. The background flux 
was subtracted from the flux in the source aperture to give 
the source flux. 



3.2 ISO-PHOT Observations 

The data were reduced using PHOT Interactive Analysis 
(PIA) Software (version 9.1)[| from the edited raw data prod- 
ucts. The sources were observed using a chopper throw of 
three arcminutes in a rectangular configuration. The ISO- 
PHOT integrations were readout and stored in a series of 
ramps containing information in time slices. Destructive 
readouts were discarded and the ramps were corrected for 
non-linear detector response. 

Measurements using ISO-PHOT could be taken in three 
main modes: staring, where the total exposure time is spent 
upon the source; maps, where the source under study and 
its immediate surrounding area are mapped by overlapping 
pointings; chopped, where the exposure time is divided be- 
tween an on-source position and an off-source (sky back- 
ground) position. In this paper we present data from the 
staring and chopped observation modes which require dif- 
ferent reduction methods. 

The staring ISO-PHOT-P observations (i.e. 25/im) were 
reduced using the standard procedure. Glitches and their 
after-effects caused by cosmic ray hits were removed. The 
ramps were then reduced to a single data point per ramp 
with further discarding of destructive readouts and glitches. 
The signal-per-ramp data were combined to give the signal- 
per-pixel. The signal was calibrated using the internal fine 
calibration source (FCS) measurement. 

Since the IR background at 25/im is dominated by zo- 
diacal light, we estimated the IR background for eac h of our 
sources using the zodia cal emission model given in 



Rowan- 



Robinson et al. (1991) evaluated at the appropriate solar 
elongation angle at the time of observation. This calculated 
background was then subtracted from the staring measure- 
ment giving the on-source flux. However, due to uncertain- 



data of sufficient resolution are not available to accurately 
subtract 25/im backgrounds for these sources. For this rea- 
son the 25/im fiux has not been included in the subsequent 
analysis but is considered an upper limit. 

Chopping is a widely used observational technique to 
observe faint sources since higher sampling of the data, from 
repeated measurements at both on and off-source positions, 
provides lower noise levels than those obtained from single 
exposures. In the case of ISO-PHOT, the chopping tech- 
nique was used to reduce the low frequency detector noise. 
The technique worked well in eliminating the long-term sig- 
nal drift of the detector since the drift timescale is longer 
than that of a chopper cycle (Abraham et al. 2001). How- 



ever, during the mission it was discovered that the chop- 
ping technique did not perform according to specifications 
due to the introduction of additional instrumental effects. 
An extensive description of the problems associated with 
this observation temp late (and their corrections) is given in 
Abraham et al. (2001). These effects include short term de- 
tector transients which, for the ISO-PHOT-C detectors, re- 
sulted in signal loss (i.e. underestimation of the source flux) 
due to the on-source integration time being shorter than the 
detector stabilisation time|^ Thus observations taken with 
the highest c hopper frequencies ar e most strongly affected 
(Figure 2b in Abraham et al. 2001). A chopper offset effect 
was also identified which can mimic or hide a real detection 
caused by inhomogeneous illumination of the focal plane. 
New deglitching routines were devised specifically for this 
observation mode since the standard algorithms, devised for 
staring measurements, were inefficient as the signal was of- 
ten not stabilised within a chopper plateau. 

The severity of the problems associated with the 
chopped mode was not known until well into the ISO mission 
by which time the PHT22 data presented here had already 
been taken|^. However, since the problems were identified 
during the lifetime of ISO, it was possible to obtain calibra- 
tion measurements in order to determine corrections. The 
analysis and c onclusions derived fro m such measurements 
are detailed in Abraham et al. (2001). They report that the 



signal loss due to chopping is dependent upon the on/off- 
source integration times as well as the flux of the source and 
background. This effect is found to be quite strong in the 
ClOO detector, where the source flux is signiflcantly under- 
estimated, and less so for the C200 detector (~ 10 — 20%). 
Abraham et al. (2001 ) describe the calibration issues and 
have derived corrections which have been implemented in 
PIA v9.1. In addition they also devised a new reduction 
method ('ramp pattern') for chopped measurements which 
is a more stable reduction method and is less sensitive to 
glitches. After the correction for the long term signal drift 



^ This aperture was chosen following a 'curve-of-growth' analysis. 
^ PIA is a joint development by the ESA Astrophysics division 
and the ISO-PHOT Consortium led by MPIA, Heidelberg. Con- 
tributing /SO-PffOT Consortium institutes are DIAS, RAL, AIP, 
MPIK and MPIA. 



* the response time taken by the detector to measure the true 
flux of a source 

^ Some 7000 measurements were taken in chopped mode i.e . a 
substantial fraction of ISOPHOT's Legacy ( Klaas et al. 200l[ ) 
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and application of new deglitching algorithms, the repeated 
on and off-source data are reduced into a single, high signal- 
to-noise 'pattern' consisting of only eight points (four on- 
source and four off-source). 

The 'ramp pattern' method was therefore used to reduce 
the data. Following power calibration the background mea- 
surement was subtracted from the on-source and the result 
was flux calibrated using the internal FCSl measurementp^ 
For the ClOO detector, it is recommended to derive the flux 
using only the central pixel (pixel 5). This is because the 
signal loss corrections have been determined solely on data 
from pixel flve (to date). 

The point spread function was larger than both of the 
ClOO and C200 detector arrays. Corrections have been esti- 
mated for the signal loss for the entire array and for pixel 

Klaas priv. 



2001 



5 of the ClOO array only ( jSalama 200C 
comm.). The appropriate factors were applied to the data 
to correct for the unmeasured flux in the wings of the PSF 
lying outside the detector array or pixel. 



3.3 Error Budget 

A signal-to- instrumental noise (comprising the statistical er- 
rors propagated through the reduction software) ratio of 
greater than three was used to determine the signiflcance 



of source detection. In addition, the uncertainty associated 
with photometric calibration must also be considered to pro- 
vide a representative total error on any measurement. The 
calibration of ISO data has been known to be problem- 
atic, but has been intensive ly investigated for all instruments 
(Metcalfe and Kessler 2001) and is now well understood. We 
consider the calibration status presented in the documents 



Klaas et al. (2001) for ISO-PHOT and Cesarsky fc Blom 



maert (2000) for ISO-CAM to derive calibration accuracies 



which are appropriate for our measurements. 



3.3.1 ISO-PHOT 

The essence of the calibration error for ISO-PHO T lies not 
only in the chopped mode strategy but also in inaccuracies 
of the measurement internal FCS calibration source. Many 
calibration tests have been carried out to quantify the lev- 
els of this error for the ISO- PHOT instruments. We use the 
photometric accuracies from Klaas et al. (200ll ) p^ . For the 
staring observations at 25^im the photometric accuracy is 
20%. However, investigations into the chopped mode pho- 
tometry show that, at faint flux levels (ClOO < IJy and 
C200 < 2 Jj/) , a careful consideration of the cirrus confusion 
noise is mandatory to determine the overall photometric ac- 



curacy o f the measurement (4.braham et al. 2001 
al. 2001|). 



Klaas et 



The data for this program was taken in a rectangular 



The FCSl measurement is thought to be more reliable than 
the FCS2 measurement, so only the FCSl measurements are used 
for flux calibration (Abraham priv. comm). 

The accuracies given in this reference are for the offline 
pipeline processing software (version 10) and are also appropriate 
to those resulting from an interactive reduction using PI A v9.1 
(Klaas priv. comm.). 



chopping mode with a single off-source position. An ISO ob- 
server could define the magnitude of the chopper throw, but 
the actual location of the background measurement was de- 
termined by the telescope configuration (i.e. roll angle) at 
the time of observation. Thus a consideration of the fiuc- 
tuations in the IR background due to cirrus or unresolved 
galaxies must be included. 

Therefore in the error budget we must consider three 
sources of error: 

• instrumental the statistical error propagated through the 
reduction software 

• photometric the error associated with the accuracy of the 
calibration measurements 

• cirrus confusion the error associated with the IR back- 
ground measured by the off source position. 
Measurements are generally instrumental noise limited but 
some may be limited by the cirrus confusion noise. 

The conf usion noise was d etermined using the prescrip- 
tion given in Kiss et al. (2001) (the equation numbers cited 



in this paragraph are references to the equations in 



Kiss 



et al.). Firstly, the zodiacal light contribution was sub- 



tracted using estimates from th e COBE/DIRBE zodiacal 
light model ( Kelsall et al. 199^ interpolated to the ISO- 
PHOT wavelengths (Abraham and Klaas priv. comm.). The 
cirrus confusion noise was then calculated using Eq uation 4 
for 60 an d 180/im and Equation 6 for 90pm from 



Kiss et 



al. (2001). The noise is also dependent upon the size of the 



chopper throw (the confusion noise increases with increas- 
ing separation between the on and off source positions) and 
thus the noise is corrected using Equation 7. Finally, a 20% 
correction is included to account for having only a single 
off-source position. 

For all measurements, bar one, the cirrus confusion 
noise is much lower than the instrumental noise i.e. the mea- 
surements are instrumental noise limited. For these mea- 
surements, the sta ndard photometric accuracies given by 
Klaas et al. (200l|) of 30% for ISO-PHOT-ClOO and 40% 



for ISO-PHOT-C200 apply (Klaas priv. comm.). For the 
cirrus-confusion noise limited measurement (TXS0052-(-471 
at 180/im) we apply a calibration error of 50% based upon 
the level of the cirrus confusion noise (Klaas priv. comm.). 



3.3.2 ISO-CAM 



The photometric ca libration accuracy is given by Cesarsky 
& Blommaert (2000). For the low fiux range of our sources, 
an accuracy of 30% is appropriate. 



4 MODELS FOR THE INFRARED EMISSION 
OF GALAXIES 

In this section we give a brief summary of the development 
of IR models, including a description of the models used to 
fit the ISO-IRAS SEDs obtained here. 



4.1 Background 

In general, the emission of IRAS g alaxies was modelled 



using a mix ture of components (e.g. de Jong et al. 1984 
Helou 1986|): 
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(a) an interstellar dust component which is heated by am- 
bient starlight associated with an old quiescent population. 
This component was named 'cirrus' since it has the same 
characteristics as cirrus seen in the Galaxy i.e. peak ing at 



lOO/^m with some emission at 60pim ( Low et al. 1984 ) 
(b) a 'starburst' component representing warmer emission 
which is Eissociated with emission from dust clouds within 
which active star formation is occurring. It d ominates the 
IR emission at 50-100/jm peaking at ~60/im (Soifer et al 
1984 Nowan-Robinson et al. 198^. 



Rowan-Robinson & Crawford (1989, hereafter RRC89) 
incorporated two improvements in these component models 
based upon information derived from IRAS colour-colour 
diagrams and source SEDs. Firstly, they used an exact solu- 
tion to the radiative transfer problem through star forming 
regions to improve the model of the starburst component. 
Secondly, they introduced a third 'Seyfert' component. This 
component contributed to the MIR emission and displayed 
a broad peak at ~ 25Mm , reflecting the characte ristic IR 
emission seen in AGN (e.g. [Edelson fc Malkan 198^ ). RRC89 
found that their IR starburst model fit the observed SED of 
the archetypal ULIG Arp 220 extremely well. 

Models comprising some combination of the compo- 
nents detailed above provide good fits to the observed SEDs 
of most IR galaxies. The contribution of each component 
to the total power varies from source to source. For exam- 
ple, some ULIGs have been observed to have IRAS colours 
similar to those of Seyferts [S25/S60 > 0.2 as deflned in |d^ 
Grijp et al. (1985)] and were defined to be 'warm' ULIGs. 
The best-fitting models to the SEDs of such sources would 
include a stronger Seyfert component than is used for mod- 
elling 'colder' IRAS galaxies. 



4-1.1 Detailed Pure Starburst Models 



Kowan-Kobmson & Jjjtstatlnou (iyya, hereafter RRE93) 
presented an improved radi ative transfer starburst model 
using the dust grain model of Rowan- Robinson (1992). They 



found their models fitted the emission of ULIGs well but re- 
quired a higher optical depth than was needed to explain 
the emission of less luminous infrared galaxies (i.e. from 

Tuv = 200 to Tnv = 500). 

'State-of-the-art' pure starburst models have been de- 
veloped by ERSOO again using radiative transfer theory. 
These models include a simple model for the evolution of HII 
regions surrounding star s, employ an advanced dust grain 
model containing PAHs ( Biebenmorgen & Kriigel 1992) and 



an evolving st ellar population derived from th e stellar syn- 
thesis codes of Bruzual & Chariot (1993, 1995). The models 
assume a fixed initial star formation rate which exponen- 
tially decays with an e-folding time of 20Myrs. Other vari- 
able parameters are the starburst lifetime and initial visual 
optical depth. The latter parameter (which has four discrete 
values Tv = 50, 100, 150, 200) also determines the star forma- 
tion efficiency of the molecular clouds which is set to be 25% 
for Tv = 50 and halves for each step in r„. The SEDs of star- 
burst galaxies M82 and NGC 6090 were demonstrated to be 
well fit by one or a combination of two of the models from 
this set. 



4.1.2 Detailed A GN Models 

RRE93 modified the MIR Seyfert model of RRC89 to model 
the emission of an AGN through an optically thick, spher- 
ically symmetric dust shell with density n{r) oc r~^. This 
model was found to fit the MIR emission of some IRAS 
galaxies well. Such Seyfert models predict the presence of 
an unobserved lO^m emission feature. RR95 constructed a 
geometric model of the obscuring material which suppressed 
this unobserved feature. He assumed the obscuring matter 
was an ensemble of isolated, approximately spherical, dust 
shells. The distribution of the obscuring material in this 
manner caused the emission and absorption of photons at 
lO/xm to be approximately equivalent since the shells act as 
good approximations to a black-body. The shells were lo- 
cated at distances from 1 to 300pc from the central source 
with density distribution n{r) oc r~^ and had dust temper- 

atures in the range of 160-1600K. 

Efstathiou & Rowan-Kobinson (1995) modelled the IR 
emission of AGN using accurate solutions of the axially sym- 
m etric radiative-transfer problem [an ad aptation of the code 
of Efsta thiou fc Rowan- Robinso n (1990)]. A multi-grain dust 
model ( Rowan- Robinson 1992 ) was used, a range of den- 
sity distributions were explored and three tori geometries 
were considered; fiared disks, tapered disks and anisotropic 
spheres. On comparison with observational data, they found 
the best-fitting models to be those with tapered disk geome- 
tries with an opening angle ~ 45° consisting of dust grains 
following density distribution of n(r) oc r~^ . Flared disks, 
anisotropic spheres and tapered disks with different param- 
eters predicted spectra with either too narrow an infrared 
continuum or strong 10fj,m features. The geometry and den- 
sity distribution of the preferred tapered disk model sup- 
pressed this 10/xm feature. Based upon these tapered disc 
models, EHY95 successfully modelled the emission of Seyfert 
galaxy NGC1068. However, an additional component of op- 
tically thin dust {Ay ~ 0.1 — 0.5), situated between the 
narrow and broad line regions, was required to explain all 
of the infrared emission. 

Alternatively, models of warped disks around quasars 
have been advocated to explain both the MIR and F IR emis- 

1989| ) . In such 



Sanders et al. 



sion of quasars and ULIGs (e.g. 

models, the NIR-MIR emission originates from hot dust clos- 
est to the central engine which reprocesses nuclear light. The 
FIR-mm emission emanates from cooler dust located in the 
outer regions of the warped disk which is also illuminated by 
the central power source. Warped disk models using radia- 
tive transfer theory are yet to be developed. RR95 demon- 
strated that the existing warped disk models were unable 
to explain all the IR emission of some PG quasars. Their 
SEDs were found to be better modelled using a sum of 'disk', 
'Seyfert' and 'starburst' components. 



4.2 Combined Models 



Green & Rowan- Robinson (1996) found that the IR emis- 
sion of the HyLIG IRAS P09104-I-4109 was weU explained 
using a flared-disk Seyfert m odel without any contribution 
from a starburst. Moreover, Glranato, Danese fc Frances- 



chini (1996) found that the emission of IRAS F10214, IRAS 



P09104 and IRAS F15307-F3252 could aU be modelled us- 
ing a heavily dust enshrouded quasar with a thick dust torus 
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model ( Sranato fc Danese 1994 ) alone. In contrast, Green 
& Row m- Robinson (199(:) successfully modelled the emis- 
sion of IRAS F10214-1-4710, using a radiative transfer model 
requiring both starburst and flared disk Seyfert components. 
RR2000 found that a combination of both starburst and 
AGN models was required to explain the SEDs of over 70% 
of all HyLIGs known at the time. 



In summary, the infrared to sub-millimetre contin- 
uum of HyLIGs may be modelled by employing elaborate 
torus/disk models for AGN. However such models do not 
conclusively demonstrate the power source of the emission 
since they often fail to explain all of the emission beyond 
50/im for which a starburst component is required. 

Of the models described above we use the following (and 
combinations thereof) to compare to our ISO-IRAS SEDs 
presented in the subsequent section: 

(1) RR95 MIR emission model for quasars 

(2) EHY95 inclined dust torus models^ 

(3) ERSOO pure starburst modelQ 



the contributions of both components to freely vary. It is im- 
portant to note that by using a combination of the models 
described in the previous section we are by no means sam- 
pling the full parameter space of starburst or AGN models. 
The number of free parameters we can use is restricted by 
the sampling of the SED. This enables us to determine the 
quality of the combined fit by calculating the reduced chi 
squared (Xr) estimator. The combined models with Xr ly- 
ing between the minimum chi squared value (Xflmin) Eind 
l+xkmin were taken to be the best-fitting combinations. 
The degeneracy within this range of combined models was 
investigated to determine the best-fitting model. The contri- 
bution levels of the starburst and AGN to the total infrared 
power are then directly obtained from the best-fitting model. 

Additionally, optical and/or NIR spectroscopic classifi- 
cation was used to verify consistency with the inclination of 
the equatorial plane of the dust torus to the line of sight. 
(The models assume a half opening angle (6') of 9 = 60°. 
Thus to obtain a narrow-line or Seyfert 2 like spectrum the 
broad line region must be obscured i.e. the inclination angle 
i must be -30° < i < 30°. For |i| > 30° QSO or Seyfert 1 
spectral signatures would be seen). 



5 RESULTS 

The fluxes and associated errors obtained per source are pre- 
sented in Table For both ISO-CAM and ISO-PHOT ob- 
servations, detections were considered to be real if they were 
at least three sigma detections otherwise a three sigma up- 
per limit is given. /SO- CAM images and SEDs are presented 
below. Where available, radio fiux contours taken from the 



FIRST survey ( [Becker, White fc Hefland 1995| ) are over- 
plotted on the ISO-CAM images , otherwise data from the 
NVSS survey (ICondon et al. 1998|) are used. The IRAS Fa.mt 
Source Catalogue (Moshir et al. 1992) one-sigma positional 
error elhpses are also shown. The ISO-CAM a.nd ISO-PHOT 
fluxes were combined with the IRAS d&ta, and any previously 
published IR, sub-mm or optical data to form the SEDs. 

These SEDs have been compared to three models of 
starbursts and Seyferts described in Section 4.2, In the first 



instance, best-fitting models from the FIR to the MIR for 
each of these three models were selected using the mini- 
mum reduced chi squared estimator. The best-fit of each 
model type outlined above are overplotted on the SEDs for 
comparison. From these plots we obtain an idea of which 
model type is the most dominant contributor (i.e. whether 
the source is starburst-like or AGN-like). Additionally, this 
plot demonstrates whether single component models alone 
are able to explain the emission across the entire IR range. 

We also combined the inclined dust torus models of 
EHY95 with the pure starburst models of ERSOO allowing 



EHY95 derive models to explain the emission of NGC 1068. 
We do not employ the NGC1068 models themselves, but we use 
the inclined dust torus models which are the basis of the NGC1068 
models (i.e. NGC1068 models without the additional conical dust 
component). 

It is important to note that the starburst models were de- 
veloped to explain the IR emission of far less extreme sources 
(L 10"'^''~^"'^Lq). Thus parameters such as star formation effi- 
ciency may be underestimated. 



5.0.1 Key 

In all cases the ISO-PHOT data is given by • and the ISO- 
CAM by ■. IRAS good or moderate quality data are given 
by Upper limits are plotted as the unfilled symbol o, □ 
and A with a downwards arrow plotted directly below. 

Photometric data obtained from literature are also in- 
cluded in the fitting algorithm. They are plotted on the 
SEDs as unfilled diamonds (with arrows below if upper lim- 
its). The source of any additional data is detailed in the 
figure captions. 

The starburst model (ERSOO) is denoted by a solid line, 
the spherically symmetric quasar and Seyfert model (RR95) 
by a long dashed line and finally the inclined torus models 
(EHY95) by a dotted-dashed line in the bottom left part of 
Figures 

In the combined model SEDs the starburst model is 
plotted as a long dashed line and the dust torus model as 
a dashed-dotted line. A solid line represents the best-fitting 
combined model. 



5.1 IRAS F00235+1024 

Fluxes are given in Table |l| and the SED is plotted in Figure 
|l[ The IRAS 60^im detection is confirmed by ISO-PHOT. 
We also obtain significant detections at 90 and 180/im. The 
comparison of the IR emission to the single component mod- 
els reveals that the emission over the entire wavelength range 
is generally more consistent with a starburst but it does not 
explain all the data well. The predicted FIR emission by 
both AGN models is substantially below that of the mea- 
sured values. Therefore none of the pure starburst or pure 
AGN models alone can fit the IR emission of this source 
well. 

The selected best-fitting combined model includes a 
starburst model of age 57Myr and with initial UV opti- 
cal depth of 880, contributing at a level of 63% to the 
IR luminosity. The remaining 37% is attributed to a torus 
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Name 


IRAS F00235+1024 


TXS0052+471 


IRAS F14218+3845 


IRAS F15307+3252 


IRAS F23569-0341 


^Va. V e 1 e n g t h 


Flux 


Flux 


Flux 


Flux 


Flux 














IRAS 12 


< 173 


< 78.03 


< 96.90 


< 45.0 


< 87.0 


IRAS 25 


< 193 


< 59.14 


< 74.90 


80.0 ± 24.0 


< 142 


IRAS 60 


428 ± 55.6 


223 ± 53.6 


< 565.0 


280 ± 27.0 


347 ± 30 


IRAS 100 


< 938 


< 1732 


< 2100 


510 ±62.0 


< 792 


CAM 6.7 


0.917 ±0.384(3.42) 


< 1.608 


0.785 ± 0.258(5.46) 


3.47 ± 1.06(17.33) 


< 0.73 


CAM 15 


6.752 ± 2.143(9.66) 


2.139 ± 0.828(4.09) 


3.228 ± 1.036(6.44) 


20.58 ±6.21(28.98) 


< 3.8 


PHOT 25 


< 376.7 


< 249.6 


< 263.9 


< 318.5 


< 963 


PHOT 60 


303.4 ± 100.3(7.19) 


< 356.5 


< 168.5 


405.9 ± 139.3(6.00) 


< 303 


PHOT 90 


477.9 ± 147.9(13.13) 


166.5 ± 73.6(3.08) 


163.5 ±61.02(4.5) 


368.4 ± 115.7(10.78) 


< 142 


PHOT 180 


804.7 ± 393.3(3.56) 


265.4 ± 158.7(5.66) 


< 1765 


413.6 ± 176.3(6.79) 


< 610 



Table 1. A compilation of IRAS Faint Source Catalog and ISO fluxes for the /50-HyLIGs. Three sigma upper limits are quoted for non 
detections. For TXS0052+471 fluxes are from IRAS Faint Source Reject Catalog not the Faint Source Catalog. Errors shown include 
calibration errors, signal to noise is given in brackets. 



model inclined at 17° to the line of sight. Models with 
XHmin<Xfl< l+Xflmin have a distribution of contributions 
with a mean starburst to AGN ratio of 64% : 36%. The 
plotted model reflects this ratio. The mode of the contri- 
bution distribution is 85 - 90% : 15 - 10% (not plotted). 
Such combinations can explain most of the IR emission but 
greatly overpredict the 90^m emission and require an AGN 
to power all of the optical emission [this is inco nsistent with 



evid ence from high resolution I-band imaging (Farrah et al 



2002 |)]. [ The distribution of AGN fraction in the low Xr range 
extends to a maximum 75% AGN contribution. Such com- 
posite models over-estimate the emission measured by ISO- 
CAM and predict a 25/im power [in logioii^S,^)] of ~ 12.3. 
Without a 20-30/im constraint we cannot categorically rule 
out a stronger AGN in this source. Nevertheless, given the 
data available, the mean of the combination distribution and 
the overpredictions found using the extreme AGN contribu- 
tions (of 10% and 75%) we conclude that it is likely this 
source is predominantly (~ 60%) starburst fuelled and an 
AGN component is required to explain full the SED. 

The dominance of the starburst contribution to the IR 
emission of this source determined from the fit of the SED 
is consi stent with observatio ns at other wavelengths. For ex- 
ample, Wilman et al. (1995) find the soft X-ray upper limits 
from ROSAT-HRI observations are consistent with an origin 
in a starburst. In addition, the predominant starburst fu- 
elling is consistent with the narrow-line classification based 
upon the optical spectrum of this object (McMahon et al., in 
prep.). Moreover, the inclination of the torus model with re- 
spect to the line of sight (17°) also supports the narrow line 
spectral classification and the non-detection of X-rays since 
both the broad line region and soft X-ray emission would be 
obscured by the dust torus in the line of sight 



The ISO, IRAS, optical [DSS and HST ( |Farrah et al 
2002|)n.nd radio (NVSS) source positions are in good agree- 
ment. 



5.2 TXS0052+471 

Fluxes are given in Table |l| and the SED is plotted in Figure 
H The IRAS, ISO-CAM and radio positions are coincident 
(lower left hand corner of the ISO-CAM image). The DSS 



image of the corresponding region reveals no optical source. 
Aperture photometry of the ISO-CAM image was therefore 
performed at this radio peak rather than the pointing centre. 
This radio position does not lie on the central pixel of the 
ISO-PHOT ClOO detector but lies between pixels 2 and 3 of 
the array. Therefore the fluxes quoted for the ClOO array are 
calculated from the entire array with the appropriate PSF 
correction applied. The 60 and 90/im fluxes therefore may 
be underestimated due to signal loss outside the array and 
in the gaps between the detector pixels 2 and 3. 

Nevertheless, these observations provide conclusive ev- 
idence that this IRAS Reject Catalogue object is a real 
infrared source. The SEDs clearly display this source has 
IR emission more consistent with an inclined torus model 
rather than any starburst model. The inclined torus model 
also satisfies the sub- millimetre limits from Ivison (RR2000, 
priv. comm.). Investigations of the combined model degen- 
eracy reveals that all low models favour a torus model 
with orientation of 26° with a — 4% starburst contribution. 
However the combined models, which include a low level 
starburst contribution, are forced to pass extremely close 
to the sub-millimetre limits and therefore a null starburst 
contribution is preferred. 



5.3 IRAS F14218+3845 

Fluxes are given in Table ^ and the SED is plotted in Fig- 
ure ^. This source was observed twice by ISO-PHOT and 
/SO- CAM under two separate programs. We have therefore 
combined the data from both observations of this source. 
In total, we have measurements at 60, 65 fim and repeated 
measurements at 6.7, 15, 90 and ISO/im. The fluxes obtained 
from both observations are mutually consistent. 

Both /SO- CAM images show a weak source close to the 
centre of the image which associated with the source seen in 
the DSS image. The /SO- CAM position lies within the cen- 
tral pixel (pixel 5) of the ISO-PHOT ClOO array. Using the 
pixel 5 fluxes we find that the IRAS 60 and 100/xm fluxes 
of 0.565 and 2.1Jy are not confirmed by the measurements 
taken with ISO-PHOT at 60 and 90/im. There is evidence 
from both the IRAS (the source is flagged as extended at 
60 and 100 ^m) and the ClOO observations (flux is recorded 



© 2000 RAS, MNRAS 001, 0-|l| 



Observations of HyLIGs with ISO 9 




IRAS F00a35+10a4 



IRAS F00a35+10a4 



enCMyr TTrT=HSO EB500 Sterturet Mod^sl 
RR9& tjuHBiM- Mod^l 
AKE InoliiieiJ Torue fl=5^+5° Mod^sl 




&7MJT TOT=BflO EB500 Sterturet Mod^sl 
AWE Inclined Torus 0=17' Mod^?l 
Totel 



12 13 

00 S6 07.07 +10 41 31,4 




12 13 

00 SO 07.07 +10 41 31,4 



Figure 1. Images for IRAS F00235+1024. Top Left /50-Cj4M image at 15/^m. The circled white pixels are the detected source emission. 
The IRAS one sigma error ellipse is overplotted (dashed ellipse) as are NVSS contours within this region. Top Right The optical DSS 
image zoomed in on the source location. The ellipse marks the IRAS error ellipse and the circle a rad ius of 22" around th e ISO source. 



Bottom The SEDs contain a combination of IRAS and ISO data. Extra I band data is taken from Farrah et al. (2002). Bottom Left 
Overplotted on the SED are the best-fitting individual models: starburst model (solid line) from ERSOO, dust torus model from EHY95 
and the MIR quasar and Seyfert model (heavy dashed line) from RR,95. Bottom Right This SED displays the best-fitting combined model 
of starburst and AGN components. The starburst component (ERSOO) is plotted with a long dashed line and the inclined dust torus 
component (EHY95) with a dashed-dotted line. The sum of the two components is plotted shown as a solid line. 



in pixels other than pixel 5) that the IRAS fluxes may be 
contaminated by cirrus. Any starburst or AGN torus could 
not be greater than one arcminute in extent, as implied by 
the unresolved IRAS size, which corroborates that the IRAS 
fluxes are contaminated by cirrus. Such contamination thus 
explains the lack of confirmation of the high IRAS fluxes 
by ISO-PHOT. However, at this stage, we cannot exclude 
the possibility that the emission detected by IRAS origi- 
nates from more than the central /S'0-Cj4Msource. The area 
around IRAS F14218-I-3845 needs to be mapped in greater 
resolution by future FIR instruments for confirmation. 

We find that at 60 and 65/im the source is not de- 
tected in pixel 5. Significant detections are obtained at 90 
and 180/im and the fiuxes obtained from the two repeat ob- 
servations are consistent within the total errors. At QO/im 
the flux given is derived only from pixel 5 as recommended 
for the ClOO array. The ISO-CAM detection and HST op- 



tical position (Farrah et al. 2002) are in agreement and lie 



within pixel 5 of the ClOQ array. We are therefore confident 
that the flux obtained emanates from the identified quasar 
at z ^ 1.21. Unfortunately at ISO/xm, the /SO- C^M source 
lies directly in the centre of the C200 array. Since the emis- 
sion from the central quasar cannot be separated from any 
cirrus (or possible additional source) contributions within 
the C200 array, we use the strongly detected 18Q/im fiux 
only as an upper limit to the 18Q/.im emission of the QSO. 

Figure ^ shows the SED obtained for the identified 
QSO. We used a wide range of torus models (varying in torus 
opening angle, line of sight orientation and normalisation) 
to test if the source emission could be explained by a torus 
model alone. In all cases good fits to the MIR emission failed 
to explain the 90/im emission and models passing through 
the 90/im detection greatly overpredicted the MIR emission 
by 2-3 orders of magnitude. Therefore a combined model was 
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Figure 2. Images for TXS0052+471. Top Left ISO-CAM image at 15/im. Tlie circled wiiite pixels are the assumed source detection. 
The IRAS one sigma error ellip se is overplotted (das hed ellipse) as are NVSS contours within this region. Positions of the peaks of the 
double radio source detected by Blundell et al. (1994 ) are plotted as crosses. Top Right The optical DSS image zoomed in on the source 
location. The ellipse marks the IRAS error ellipse and the circle a radius of 22" around the ISO source. Bottom The SEDs contain 
a combination of IRAS and ISO data. Extra sub-millimetre limits were taken from RR2000 (reference therein Ivison priv. comm.). 
Bottom Left Overplotted on the SED are the best-fitting individual models: starburst model (solid line) from ERSOO, dust torus model 
from EHY95 (dashed-dotted line) and the MIR quasar and Seyfert model (heavy dashed line) from RR95. Bottom Right This SED was 
generated as the best-fitting from the combined model analysis. This shows that the best-fitting was provided by solely an AGN torus 
with no starburst contribution (solid line). 



preferr(d. As above, the components are allowed to freely 
vary in contributions. The best-fitting models (i.e. those be- 
tween XRmin<XR< l+Xflmin) displayed a common feature 
where the MIR-optical emission is dominated by the torus 
emission and the 90/im flux explained by a starburst. De- 
generacy in the low Xr range {xRmin<XR< was 
investigated revealing over 80% of the combined models pre- 
ferring a power ratio of 75 — 85% : 15 — 25% starburst:torus. 
The best-fitting model shown in Figure ^ reflects this pref- 
erence. The combined model consists of a 1.7Myr starburst 
with an UV optical depth of 880 contributing at a level of 
74%, the remainder being accounted for by an torus model 
inclined at 59.5° to the line of sight. This orientation of the 
torus model with opening angle of ~ 120° implies the broad 
line region is visible. The model is therefore consiste nt with 
the detection of broad lines in the optical spectrum (Farrah 



et al. 2002). The best- fitting models also display a slight 



over-prediction of the 15^m emission. It is most likely this 
is due to the fact that we are not sampling the full parameter 
space of starburst and AGN models. 

Within the low Xr range we can find solutions with 
AGN contributions as low as 10% and as high as 40%, but for 
both of these extremes the combined model overpredicts the 
MIR emission detected by ISO-CAM. The conclusion that 
the thermal emission is predominantly of starburst origin 
remains unchanged. 

Despite the ISO flux being almost a factor of 10 lower 
than those detected by IRAS the ISO detected fluxes and 
fits confirm that the QSO (at a redshift of 1.21) remains a 
bona fide HyLIG. 
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Figure 3. Images for IRAS F14218+3845. Top Left Co-added ISO-CAM image at 6.75Aim. The circled white pixels are the detected 
source emission. The IRAS one sigma error ellipse is overplotted (dashed ellipse) as are NVSS contours within this region. Top Right 
The optical DSS image zoomed in on the source location. The ellipse marks the IRAS error ellipse and the circle a r adius of 22" around 
the ISO source. Bottom The SEDs contain a combination of IRAS and ISO data. Extra I band data was taken from Farrah ct al. (2002) 
and R, and B magnitudes from the APM survey. Bottom Left Overplotted on this SED are the best-fitting individual models: starburst 
model (solid line) from ERSOO, dust torus model from EHY95 (dashed-dotted line) and the MIR quasar and Seyfert model (heavy dashed 
line) from RR95. Bottom Right This SED displays the best-fitting combined model of starburst and AGN components. The starburst 
component (ERSOO) is plotted with a long dashed line and the inclined dust torus component (EHY95) with a dashed-dotted line. The 
sum of the two components is plotted in a solid line. 



5.4 IRAS F15307+3252 



after the formation of the star when unabsorbed Ught begins 
to leak-out both vary. 



Fluxes are given in Table |l| and the SED is plotted in Figure 
^ This object has strong ISO-CAM detections. The IRAS 
and ISO fluxes are in good agree ment. Since this s ource was 
observed using ISO-CAM-CYF ( |Aussel et al. 1998| ) we could 
use this data to further constrain the models. This additional 
data also enabled the use of a wider range of starburst (Ef- 
stathiou priv. comm.) as the number of data points greatly 
exceeds the degrees of freedom. This extended set of star- 
burst models varies not only in starburst age and r„ but also 
includes a varying time constant for the exponentially decay- 
ing star formation rate. In addition, the fraction of optical 
light that is allowed to 'leak' out from the giant molecular 
'cloud plus star' systems without absorption and the time 



Using this extended set of starburst models together 
with EHY95 dust torus models with a half-opening angle of 
6 = 30° , we obtained a moderately good fit for the emission 
of IRAS F15307-I-3252 from the optical to the millimetre. 
The strongest constraint on the fitting in the FIR is the 
1250^m limit as the high 180/xm fiux forces models to lie 
close to this limit. Investigating the parameter space be- 
tween Xfl™m<XH< l+Xflmin iudicatcs that a ~ 30% star- 
burst contribution to the IR power is favoured by the ma- 
jority of combinations. Hence we select a combined model 
where this preference is refiected and has a low Xr- 

The fitting of this source has proved difficult since the 
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Figure 4. Images for IRAS F15307+3252. Top Left ISO-CAM image at 15/^m. The circled white pixels are the detected source emission. 
The IRAS one sigma error ellipse is overplotted (dashed ellipse) as are FIRST contours within this region. Top Right The optical DSS 
image zoomed in on the source location. The ellipse marks the IRAS error ellipse and the circle a rad ius of 22" around the I SO source. 
Bottom The SEDs contain a combination of IRAS and ISO data. Extra optical data was taken from Cutri et al. (1994) and Farrah ct 
al. 



|). CAM CVF data was taken from Aussel et al. (1998) and includes the sub-millimetre upper limit from reference therein Lis 
(priv. comm.). Millimetre limits are taken from YS98. Bottom Left The SED generated combining the IRAS and ISO data together. 
Overplotted are the best-fitting model starburst model (solid line) from ERSOO, dust torus model from EHY95 (dashed-dotted line) and 
the MIR quasar and Seyfert model (heavy dashed line) from RR95. Bottom Right This SED displays the best-fitting combined model 
of starburst and AGN components. The starburst component (ERSOO) is plotted with a long dashed line and the inclined dust torus 
component (EHY95) with a dashed-dotted line. The sum of the two components is plotted with a solid line. 



non-de tection o f PAH features in the CVF spectrum (Aus- 



sel et i 



1. 1998) greatly limits any starburst contribution. 



Despite this, a strong starburst component is required to 
explain the 180/im emission since any of the torus or spher- 
ically symmetric AGN models cannot explain all the FIR 
emission beyond 60^m. Therefore in addition to the X^min 
to l+Xflmin hmits we also include that PAH features should 
not be significantly detectable above the continuum in the 
NIR. A starburst contribution is required to explain the FIR 
emiss ion which must also satisfy the sub- millimetre (taken 
from Aussel et al. 199^ and millimetre (Yun & Scovilk 



This 



is widely seen in AGN where PAHs are destroyed by the 



intense radiation field. 



1998, hereafter YS98) limits. Figure ^ shows an example 
of such a combined model, however the CAM 15/^m and a 
few of the higher wavelength ISO-CAM-CYF data are over- 
predicted by the combined model by a factor ~ 2 — 3. This 
predicted excess could be attributed to either the calibration 
of the ISO-CAM data, or to the fact that the set of starburst 
and AGN models do not sample the full parameter space. 
This excess is an unavoidable prediction of the AGN mod- 
els which are required to fit the NIR-MIR data points and 
to suppress the PAH features. Therefore all the best-fitting 
combined models for this source display this 15^m excess. 
To eliminate this problem further development of the models 
is required which is beyond the scope of this paper. 

Approximately 90% of the combined models consist of 
at least a 55% AGN contribution. The remaining 10% of 
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composite models favour combinations with only a 10-15% 
AGN contribution. However such models invoke the presence 
of strong PAH features in the MIR due to the strength of 
the starburst over AGN model. Due to the lack of detection 
of such features we exclude these combinations as justifiable 
alternative fuelling compositions. Therefore our conclusion 
that this source is predominantly fuelled by an AGN in the 
infrared remains unchanged. 

NIR spectroscopy (Evans et al. 199^ ), UV/blue spec- 
t roscopv ( Sutri ct al. 1994[ ) and optical spectropolari metry 



(Hines et al. 1995) classify this source to be a Seyfert 2. Hine; 



et al. ( _995) found 13% of the optical flux to be polarised 



and detected a broad MgH line in the polarised continuum. 
Both of these results suggest the presence of an AGN with 
its light scattered into the line of sight by matter. Also upper 
limits (7.2 x 10""erg cm"^ s"^ in the 2-10 keV band) from 
soft X-ray data cannot rule out the presence of a heavily 
obscured AGN (Ogasaka et al. 1997). This limit was low- 
ered by a factor of four in luminosity (~ 4 x W^^ erp 



in the 0.1-2.4keV band) using ROSAT HRI dat a (Fabian 
et al. 1 )9(: ). If an AGN is present in this source, Fabian et 
al. (19q|)" state that it must be obscured by matter of ex- 
tremely high column density which is X-ray absorbing and 
is Thomson thick {Nh > lO^** cm~^). How ever, the observed 



polarisation of optical light in this source ( Hines et al. 199 



implies that X-rays emanating from the AGN would also 
be scattered into the line of sight unless the medium does 
not contain any free electrons (i.e. the scattering material 
is dust) or the source is exceptionally weak in X-rays. This 
then implies that there is little ionised gas in the line of sight 
which is consistent with low gas masses detected from CO 
observations (YS98) of Mg < 5 x W^h^^MQ. Combining 
the measured millimetre flux limits with the existing IRAS 
data, YS98 use the prescription of Hildebrand (1983) and 



best-fitting black-bodies to determine a dust mass within 
the range Md = 0.4-1.5 x lO^/ifg^Mo for this source. YS98 
comment the IR-sub-millimetre emission properties of IRAS 
F15307+3252 are quite different to Arp220, commonly re- 
ferred to as the archetypal starburst ULIG, suggesting that 
a starburst explanation alone is insufficient and the presence 
of a energetically significant, obscured AGN is probable. 

Overall our results are in agreement with the presence 
of a strong AGN in this source and corroborate previously 
published evidence at other wavelengths. The best-fitting 
model has starburst-to-AGN luminosity ratio of 32%:68%. 
The best-fitting torus model is inclined to the line of sight by 
i = 50° and has a half-opening angle 6 — 30°. This orienta- 
tion implies the broad line region is obscured which is con- 
sistent with the Seyfert 2 spectroscopic classification. The 
SED of this galaxy has also been successfully mode lled using 
a simpl e dust torus model at an orientation of 45° ( Granato, 
Danese|fc Franceschini 1996). This modelling wa s performed 



upo n IR AS data and optical/NIR photometry ( Cutri et al, 
1994Qlone. We obtain a similar inclination to the line of 
sight as ICranato, Danese fc Franceschini (1996| ) but a star- 
burst contribution is definitely required to explain all of the 
FIR emission. This contributing best-fitting starburst model 
has an age of 72Myr, an initial UV optical depth of 440, a 
starburst e-folding constant of lOOMyr and leak fraction of 
1.0 beginning 39.8Myr after the formation of the stars. 

Liu, Graham & Wright (1996) propose that this galaxy 
may be a giant elliptical in the process of galactic cannibal- 



ism. However, they also highlight that aspects of the mor- 
phology are reminiscent of gravitational lensing. However, 
this view is not supported by re cent high resolution imaging 
using WFPC2 on the HST (see [Farrah et al. 2002|). 



5.5 IRAS F23569-0341 

The IRAS 60/im flux could not be confirmed and there were 
no detections in the remaining filters with both ISO-CAM 
and ISO-PHOT. Hence, upper limits for this source are pre- 
sented in Table 1 but a SED is not presented. 



5.6 Luminosities and Dust Masses 

Bolometric luminosities and dust masses have been com- 
puted from the best-fitting starburst model (parameters 
given in Table P). Luminosities were calculated using the 
standard expressions for luminosity distance {Dl, Equation 
|l|) and luminosity [Lir, Equation]^ where S\ is flux. 

^ ''o^+(qo-l)[(l + 2zgo)'/'-l]} (1) 



D 



c 

Ho qo 



Ltr = 4nDi 



looo/im 



S\d\ 



i/xm 



Dust masses have been calculated usin g Equation 
from 



(2) 



taken 



Green & Rowan- Robinson (1996) with input values 



taken directly from the best-fltting starburst model. 
47r „2 S^obs "^firrii 



M, 



3 X lQ^{l + z 



■Dl-^ Tu 



(3) 



(assuming no magnification). Dl is the luminosity distance 
calculated assuming Q = 1 and Ho = 50h5okms~^ Mpc~^ 
from Equation |l|. Svobs is the flux of the model fitted to the 
observational data and S^e,„ is the flux predicted by the 
(unfitted) model at a distance of a cloud radius from the 
central source. However, since the models are calculated at 
a fiducial distance of 100 times the cloud radius, a factor 
of 10"* is applied to the calculated dust mass. Tuv is the 
optical depth at lOOOA at four discrete values of 220, 440, 
660 and 880 (where Tuv = 4.4ri,). Summations are made 
over i grain types for a given fraction (fi), mass (m^) and 
extinction cross-section at lOOOA {Cl^t.uv) to calculate the 
effect ive density of dust grains in the cloud (values are taken 
from siebenmorgen & Kriigel 1992, and Efstathiou, priv. 



comm.). 

We estimate effective temperatures and emissivities of 
the FIR emitting dust by fitting a modified black-body spec- 
trum Equation W to the FIR tail of the starburst spectrum. 
The optimal parameters for the modified black body were 
determined by using an iterative curve fitting routine until 
the resultant chi-squared values converged giving an optimal 
set of input parameters. The freely varying parameters are 
normalisation (N), Temperature (T), emissivity parameter 
(/3) and reference wavelength (Ao) 



B(A,T) 



2hc 



A3 







-[1 



(4) 



For IRAS F00235-f 1024, IRAS F14218-^3845 and IRAS 
F15307+3252 with starburst contributions, the calculated 
dust masses are > 10* assuming no magnification due to 
lensing. The mass to luminosity ratios for all sources with 
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Name RA DEC z Age r^v i Mjust L^OT lSBT ^eff 

(J2000) {J2000) (Myr) ° (Mq) (Lq) (Lq) (Lq) 



IRAS F00235+1024 00 26 04.86 +10 42 45.3 0.58 57 880 

TXS0052+471 00 54 59.95 +47 25 59.4 1.93 

IRAS F14218+3845 14 23 58.66 +38 33 12.1 1.2 1.7 880 

IRAS F15307+3252 15 32 46.97 +32 44 03.9 0.93 72 440 



11 1.28E+09 2.29E+13 1.45E+13 8.42E+12 27.3 1.68 

26 - l.lOE+14 - l.lOE+14 

59.5 2.25E+09 2.27E+13 1.68E+13 5.94E+12 24.06 1.79 

50 1.41E+08 4.79E+13 1.54E+13 3.24E+13 30.2 1.72 



Table 2. A table of the ISO observed HyLIGs. Positions are from the IRAS FSC (or FSR in the case of TXS0052+471). t„ and starburst 
lifetimes (Age) are the parameters from the best-fitting starburst model. The procedure to obtain the dust masses and luminosities of 
the sources are explained in Section 5.6. 



starburst contributions [logio{Mdust / L'f^'^)] are -4.25, -4.00 
and -5.53 for IRAS F00235+1024, IRAS F14218+3845 and 
IRAS F15307+3252 respectively. These ratios are broadly 
cons istent with the ratio found for IRAS F10214+4724 of - 
4.6 ( |Green fc Rowan-Robinson 199e| ) . The effective FIR tem- 
peratures calculated from the modified black-body fits are 
indicative of large quantities of dust at low temperatures. In 



additio: i the values of ettective dust emissivity is consistent 
with values expecte d tor local starburst galaxies (~ z e.g. 



roborating the increase in AGN-like LIGs with increasing 
luminosity relation. 

On the basis of the idea that the number of AGN 
increases in samples o f ULIGs with increasing luminosity, 
Taniguchi et al. (1999 ) postulate that it is highly probable 
that ULIGs are simply the heavily dust enshroud ed phase 
of the formation of a quasar (as was proposed by gander: ' 



Calzetti et al. 2000) 



For IRAS F15307+3252 the calculated dust mass from 
the best-fitting starburst model is 1.41 x 10^ h'^^ Mq which 
is consistent with dust mass predicted by CO observations 



et al. 1988a| ,pD- In these scenarios the AGN begins to form 
following a merger during which time a dust enshrouded 
starburst phase dominates. After sometime the screen of 
extremely thick dust, which shields the emission from the 
forming AGN, begins to break down and the active nucleus 
becomes exposed. It is postulated that this p hase quickly 



(YS98, 0-7875 - 3.375 x 10^ hjM@). The dust masses calcu- evolv es into the optical quasars seen to-day (|Lutz et al 



lated bt YS98 were based upon black/grey-body tits to the 
IR-sub-millimetre emission (IRAS data plus sub-mm lim- 
its). The dust ma sses derived from these black-body fits 
(prescription from Hildebrand 198J) are generally not ac- 
curate since they often underestimate the MIR emission if 
the fits are normalised to the FIR emission and therefore do 
not include the mass of the dust responsible for the MIR 
emission. Also the prescription is heavily dependent upon 
the value of emissivity chosen. As YS98 indicate, the result- 
ing dust masses calculated based upon the emissivity (and 
therefore mass opacity coefficient) may be uncertain by at 
least a factor of three. Nevertheless, the dust mass calcu- 
lated from our best-fitting combined SED model lies within 
the dust mass range given by YS98 yielding a gas-to-dust 
ratio of < 140.4. 



6 DISCUSSION 



6.1 



The Starburst-AGN Controversy for ULIGs 
and HyLIGs 



Many autho r s (e.g. Armus, Heckman fc Miley 1989; [Veilleux 



et al. 1995|; ganders fc Mirabel 1996| ; ^hier et al. 1996t 



Veilleux, Kim & Sanders 1999) report a general increase in 
the number of Seyferts in ULIG samples with increasing lu- 
minosity samples. If this trend continues beyond Lfir > 
10^"^ then it may be expected that a higher fraction of 
AGN-like objects in HyLIG than ULIG samples. Moreover, 
RR2000 found that in a sample of HyLIGs selected from un- 
biased surveys approximately 50% have IR emission which 
is predominantly fuelled by an AGN. This fract ion is higher 
than that found in ULIG samples [e .g. 15-25% ( Gtenzel et al. 
199i ; [Veilleux et al. 199E| ), 20-30% ( [Lutz, Veilleux fc Genzej 



199S| ). 

Barvainis et al. (1995) comment upon the similarity be- 
tween the 20/im to sub-millimetre emission of the Clover- 
leaf quasar (H1413+117, z=2.558) and F10214 suggesting 
that they are the high redshift counterparts of narrow and 
broad line AGN's respectively difi'ering only in viewing an- 
gle. HyLIGs are proposed to be the missing Seyfert 2 ana- 
l ogues to qu asars, so called 'misdirected' quasars or 'QS0-2s' 
(Hines 1995), where a QSO is seen if the pole lies along the 
line of sight and a HyLIG at any other orientation. 

Observations of HyLIGs in the hard X-ray regime pro- 
vide an obscuration independent test of the presence of an 
AGN since they have sufficient energy to penetrate the ob- 
scur ing dust torus [ a prerequisite for unified AGN models 
(see Antonucci 1993, for a review)]. For ULIGs, pales & Ar- 



naud (1985) found Arp 220 and Mrk 231 to be underlumi- 



nous in the X-ray from Einstein data and attributed this to 
absorption by dust in the line of sight. Non-detecti ons of nine 
ULIGs in the HEAO A-1 database ( |Rieke 1988| ) at higher 
X-ray energies (i.e. less sensitive to obscuration) required 
alternative explanat ions since th e interstellar absorption of 
hard X-rays is low. Rieke (1988 ) proposed three possibili- 
ties; the presence of an X-ray underluminous AGN, delayed 
AGN activity or a starburst origin for the X-ray emission. 

Typically, weak X-ray emission is also seen in HyLIGs. 
IRAS F 102 14+4710 was we a kly detected by th e ROSAT 
PSPC ( [Lawrence et al. 1994| ). [Fabian et al. (1994[ ) did mea- 
sure significant X-ray emission from ASCA observations of 
Hy LIG IRAS P09104+4109 w hich was subsequently found 
by [Fabian fc Crawford (1995[ ) using ROSAT's HRI to orig- 
inate from a cooling flow around the galaxy rather than 
scattered fight of an embedded AGN. Fabian et al. (199f) 



199£ ) if LINERs are attributed to the HII class] thus cor- 



found that IRAS F15307+3252 was undetected in X-rays 
(using ROSAT HRI) which suggested that either an un- 
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usually small fraction of the total power of a possible em- 
bedded AGN is emitted in the X-rays or that little nuclear 
X-ray flux is scattered into the line of sight by electrons. 
IRAS F15307-I-3252 was not detected by ASCA but that 
the ratio of the X-ray upper limit to infrared luminosity 
was consistent with th e presence of a highly obscured QSO 
( lOgasaka et al. 1997| ). Moreover, for four HyLIGs; IRAS 
F00235-H024, IRAS F12514+1027, IRAS F14481-h4454 and 



IRAS F14537-I-1950, [Wilman ct al. (1998D detect no con- 
firmed X-ray emission from any of the sources. It is interest- 
ing to note that even for F12514-f 1027 and F14481+4454, 
which have optical spectral signatures consistent with a 
Seyfert 2 classification, the upper limits obtained are only 
consistent if the active nucleus is atypically weak in X-rays 
or is obscured by Nh > 10^^ cm~^. Their soft X-ray emission 
bears m ore resemblance to S eyfert 2-starburst combination 
sources (Wilman et al. 1995). 

The CO luminosity (thus molecular gas) measured 
for SMMJ02399-0136 and SMMJ14011-f 0252 two sub- 
millimetre selected hyperluminous galaxies were found to 
be consistent with a substantial fraction of the infrared 



luminosit y hav ing an origin in star formation (Frayer et 
al. 199|, |l999| ). In addition, millimetre emission has been 
measured from F10214-(-4724 (|Brown fc vanden Bout 1991 



Solomon, Downes & Radford 1992) via CO molecular lines 



On the other hand, IRAS F15307-I-3252 was undetected in 
CO and at rest-frame 650/im (YS98). The derived limits on 
the molecular gas mass and the gas-to-dust ratios are below 
those typical for gas rich infra red spirals but with in observed 
ranges (YS98). Additionally, [Evans et al. (1998[ ) found that 
HyLIGs IRAS F15307-h3252 and IRAS P09104-h4109 pre- 
sented NIR emission line ratios consistent with those ob- 
served for Seyfert 2 galaxies. From non-detection of H2 they 
determined that the upper limits (< 1 - 3 X 10^° Mq) on 
the mass of hydrogen gas to be less than most gas rich IR 
galaxies. These galaxies also presented the most extre me in- 
frared/ CO l uminosity ratios known (~ 1300 — 2000). Evan^ 
et al. (L998) presented this data, in combination with the 
warm infrared colours, to be indicative of heating of a small 
amount of dust located close to the AGN. Also, for both 
IRAS F15307-^3252 and IRAS F10214, YS98 noted the FIR 
luminosity to dust mass ratio to be four times larger than 
that measured for Arp 220 and thus concluded the values 
were larger than reasonably expected for starburst domi- 
nated sources. 

The presence of an AGN in IRAS P09104+4109 and 
IRAS F15307-I-3252, is corrobor ated by det e ctions of broad 
line regions in polarise d light (Hines 1991; Hines & Wills 
1993; H|ines et al. 1991: ). Also, the detection of a highly po- 
larised giant reflection nebula in IRAS P09104-I-4109 im- 
plies that if the object were viewed from either pole it 
would be indisti nguishable from typical luminous QSOs 
( iHines etal. 1999| ). WFPC-2 and NICMOS imaging of IRAS 
P09104+ 4109 and I RAS F15307+3252 reveal bipolar mor- 
phology (Hines 199S ), again adding credence to the idea that 
some HyLIGs contain obscured AGN where emission from 
the central engine can only be observed as light scattered 
into our line of sight. 

In this paper we show that the reality of all the ob- 
served HyLIGs, except IRAS F23569-0341, has been con- 
firmed by ISO. The total far-infrared luminosities calculated 
from the best-fitting models confirm these sources are in- 



deed members of the hyperluminous class and are amongst 
the most luminous sources known in the Universe. For all 
four detected sources our results are consistent with previ- 
ously published results at different wavelengths. With so few 
sources it is difficult to extrapolate our findings to the entire 
population. Nevertheless, we can state that ALL of our 
sources require contributions from an AGN com- 
ponent to completely explain their IR SEDs. This 
suggests that obscured AGN and the hyperluminous phe- 
nomenon are linked. For IRAS F00235-fl024 and IRAS 
F14218+3845 this contribution is below that of the star- 
burst, while for IRAS F15307+3252 and TXS0052-h471 it 
is the dominant or sole contributor respectively. The high 
fraction of AGN dominated systems in this subset of a sam- 
ple of galaxies, limited only by observational constraints, 
supports the postulates of increasing AGN-like sources with 
IR luminosity. Additionally, we can hypothesise that higher 
sampling of HyLIG SEDs will reveal the need for an AGN 
component to explain all of the IR emission (as was seen for 
all our sources). We therefore conclude from the results of 
this study, coupled with the multiwavelength data available 
to date, that it is likely that HyLIGs contain AGN which 
contribute to the infrared emission. However, they are not 
energetically dominant in all HyLIGs. 

6.2 Star Formation Rates 

Accurately quantifying the degree of obscured star forma- 
tion is crucial for accurate represen tations of the star for- 
mat io n history of the Universe (e.g. |R£wan-R£bm£on^^e^^ 



19971; iCalzetti 1997|; [Hughes et al. 1998| ; [Pettini et al. 199^ ; 
Calzetti fc Heckman 199c| ). The SFR has been calculated 



from Ha, optical, UV and FIR luminosities using several 
approaches. The theory behind the estimators is that the 
infrared emission is grain re-radiation of UV and optical 
light from the photospheres of young, massive stars. Hence 
as the number of young, massive stars increases so does the 
infrared luminosity. In fact, the massive stars are also prone 
to becoming supernovae and hence the infrared luminosity 
scales with the radio. Therefore the same is true for the star 
formation rates calculated from the radio and infrared lu- 
minosities. The contention in the evaluation of the infrared 
star formation rate lies in the fact that the contribution to 
the infrared l uminosity is not solely from the rn assive young 
stars. In fact, Lonsdale Persson & Helou (1987) argued that 



stars with masses less than five solar mass are significant 
contribu tors to the FIR emission from dis k galaxies. Addi- 
tionally, Bothun, Lonsdale, & Rice (1989) highlighted that 
old disk stars can play a significant part in the heating of 
grains in normal disk galaxies. 

Currently, the majority of estimators used are based 
upon a selected initial mass function (IMF) since the star 
formation can then be inferred from the number of stars 
greater than a given mass and the power they emit. The 
crux of the determination of SFR is that the lifetimes of 
massive stars, forming the basis of the radio and IR emis- 
sion, are much shorter than a Hubble time hence the ra- 
dio/IR luminosities a re directly pro portional to the recent 
star formation rate ( [Condon 1992 ). (This proportionality 
is thus also reflected in the linear radio/FIR correlation.) 
Hence, integrating the mass implied by the IMF between 
given mass limits gives the mass of stars contributing to the 
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luminosity for starburst time period per unit area. Normalis 
ing this relation per starburst lifetime gives the mass of stars 



pro duce d per year i.e. the star form ation rate (see Condon 
IQQs Q 'ollowing the prescription in [Rowan- Robinson et al 



(1997D, |the star formation rate (SFR) for the HyLIGs have 
been derived using the following approximations: 



,all 



= 1.5 X 10-^'\(t>/e)LBOL,FIR 



(5) 



adopting (f> = 1 for a Salpeter IMF and e = 2/3 for the 
effici ency of conversion between optical/UV photons to IR 
(e.g. [Pevereux fc Young 199o[ |Xu 199C| ; [Calzetti et al. 2000| ) . 
In this case we use the bolometric luminosity from the inte- 
gration of the starburst model therefore eliminating the in- 
accuracies associated with bolometric corrections. Addition- 
ally, the rates calculated are inversely proportional to the 
square of Hq , therefore the SFR is quoted in units of h^o ■ The 
star formation rate for TXS0052-f471 has not been calcu- 
lated since our results show that its IR power is derived from 
an AGN rather than starburst. The calculated star forma- 
tion rates are 10^ '^\ lO^ '^® and 10^-^* MQyr'^h^^ for IRAS 
F00235-H1024, IRAS F14218+3845 and IRAS F15307+3252 
respectively. 

6.3 Radio-IR Luminosity Correlation for HyLIGs 

A tight correlation exists, over several orders of magni- 
tude of IR flux, between FIR and rad io emission (e.g . 



Helou, Soifer, fc Rowan- Robinson 1985 ; Soifer, Neugebauer 



&i Houck 1987t [Condon 1992^ ganders fc Mirabel 1996D fo r 



Helou et al. (198! 



normal starburst and Seyfert galaxies 
parameterised the correlation by the logarithmic measure 
q=\og{SFiR/ Si^^radio) wMch has a mean value of ^2 .34 and 
a range of 2.0-2.6 for starbursts (Condon et al. 1991). If it is 
accepted that ULIGs have radio emission of starburst origin 
[i.e. synchrotron emission from relativistic electro ns acceler- 
ated in SNRs forming within starburst regions f Harwit & 



Pacini 975)] then it is to be expected that the FIR and ra- 



dio fluxes are correlated since they have a common source of 
emission. On the other hand, radio loud quasars and radio 
galaxies have q values between ~0-l (e.g Golombck, Miley, 
&i Neugebauer 1988^ Knapp, Bies, & van Gorkom 199^ m- 
pey & Gregorini 19931; IMarx et al. 19941; |Roy et al. 1998p. 



Whereas radio quiet quasars follow the same IR-Radi o cor- 
relation as normal star forming galaxies and ULIGs ( |Sopp 



&i Alexander 1991 



core s (iRov et al. 199^ ) 



as do Sevfert galaxies without compac t 



Smith, Lonsdale, fc Lonsdale (199^ ) 



and Sanders (1999 ) interpret this relation between RQQs 
and ULIGs to add weight to the unification of the two 
sources, but it is probably only indicative that the radio 
emission from RQQs originates from star forming regions 



rather than the active nucleus ( Bopp fc Alexander 1991 



Cram, Worth, fc Savage 1392 [~ Condon et al. (1991D also 
used the dispersiuii ill Lhe radiu/iiifrared currelaLiou Lo dis- 
tinguish between starburst (tight correlation) and AGN-like 
(more disperse) sources. For example, the RMS scatter for 
a range of galajcy types and samples displaying star forming 
activity was det ermined to be <0.2 with a mean q value of 
2.34 at 1.4GHz (Condon 1992, and references th erein) 



Rad io detections from the NVSS catalogue (Condon et 



al. 1998) were sought for the objects in the originally pro- 
posed /SO-HyLIG sample (see Table ^). An infrared source 
was regarded as being associated with a radio source if the 
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Figure 5. Radio-FIR Correlation: The figure above shows the 
correlation between the calculated far infrared luminosity be- 
tween 42. 5-122. and the radio power at 1.4GHz taken from 
the NVSS survey. Upper limits are shown for non-detections. Pre- 
viously published radio-FIR correlation are plotted including the 
sample of Markarian Galaxies (segregated by type) from Bicay et 
al. (1995, B95 ) and the samples o f luminous and ultralum inous 
galaxies from London et al. (1991, , denoted as C91) and Stan- 



ford et al. (200q , , denoted as SOO). The HyLIGs agree well with 
the ULIGs both of which have higher infrared-radio ratios than 
the Markarian Galaxies. 



one sigma IRAS and NVSS error ellipses overlapped. Figure 
|5| shows the resultant FIR-to-radio correlation. Many of the 
sources were not detected in the survey and hence are shown 
as upper limits i.e. less than 2.7mJy (corresponding to the 
completeness limit of the NVSS survey). For the HyLIGs 
the 1.41GHz flux was k-corrected in the form of (1 + z)"'^ 
where a — —d{lnS) / d{lnu) assuming a value of a = 0.7 for 
the median spectral in dex between 1.49 an d 8.44GHz (as 
calculated for ULIGs in Condon et al. 1991). 



The far-infrared luminosity from the IRAS fluxes alone 



generally calcu lated using the prescr iption given in Helou 



et al. (19851 ) and [Lonsdale et al. (198^ ) 



Lfir = 47rLiil.26 x 10""(2.58 x /so + hoo)[Wr 



(6) 



ively I 

tween 42.5- 122. 5/^rnrj. Since we have higher sampled SEDs 
and well fitting modelling the far infrared luminosities for 



A bolometric correction may be applied (C^o; ~ 1.6). How- 
ever, we use Lpjn in this Section and therefore no bolometric 
correction is applied to the result of Equation M 
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the /SO-observed sources were calculated from the inte- 
gration of the best-fitting_ model between 42. 5-122. 5/im to 
mimic Lfir in Equation pi For the remaining HyLIGs with 
well sampled SEDs, but without ISO data (see RR2000 Ta- 
bles 1-4), it was possible to perform similar SED fitting and 
then determine the far infrared luminosity as for the ISO- 
observed HyLIGs. For sources with no additional photomet- 
ric information and with a good or moderate quality IRAS 
detection at lOO/^m the far infrared luminosity was calcu- 
lated using Equation ^. However, for sources with only an 
upper limit at lOOpim the luminosity was calculated over the 
integration of a model consistent with available data and 
satisfying all the upper limits. In this manner the Lfir ob- 
tained would be more reliable than simply using Equation 
^ with a 100/im upper limit which would clearly overesti- 
mate the far infrared luminosity. The method used for each 
source is indicated in Table ^ (Lfir is only quoted for the 
KO-HyLIG sample). 

Figure |5| also shows a selection of published 1.4GHz to 
FIR luminosity ratios. The sample of Bicay et al. (1995) con- 



tains Markarian galaxies for which the spectroscopic types 
are known and includes both starbursts and Seyferts. From 
this data the starbursts display a tight correlation whereas 
the S eyferts are more dis perse. For forty LIGs and ULIGs 
from Sondon et al. (1991, also plotted) all but one have ra- 
dio emission consistent with compact starbursts. They have 
Lfir/ Pi.il ratios greater than those expected at lower lu - 
minosities. The FIRST-ULIG sample ( [Stanford et al. 200C| ) 
is als o overplotted and is in good agreement with the data 
from Condon et al. (1991). Figure H displays that the trend 
is continued for the higher redshift HyLIGs. Marked on the 
plot are the objects for which ISO data is presented and 
also the optical spectroscopic classification from various in- 
vestigations (see RR2000 Tables 1-4 and references therein). 
The two narrow line objects occupy a similar region of the 
plotted region however one has only a 1.4GHz upper limit. 
The AGNs are more disperse. 

For the sources without 750 detections or poorly sam- 
pled SEDs, the FIR luminosity calculated from the m odels 
or from the prescription given in Lonsdale et al. (1985) may 
be inaccurate due to a lack of good quality data points. In 
order to eliminate these uncertainties from the radio-IR cor- 
relation we consider only the IRAS GO/^m fiux to investigate 
the Radio-FIR correlation as most known HyLIGs sources 
have at least good/moderate detections in this bandQ. In 
addition we complement the sample with further HyLIGs 
with 60/im emission included in RR2000. The total sample 
is given in Table ^ and the resulting data is shown in Figure 
^ If the objects have ISO detections, the rest-frame 60/im 
luminosity extrapolated from the best-fitting model is used, 
otherwise a k-correction is used of the form (l-|-«)"~^ where 
a = —d{lnS) / dilnv) and is estimated to be two for starburst 
galaxies. 

It is clear from Figures I and I that the radio-IR corre- 
lation extends to higher radio and IR power than has been 



don et 



1. (1991 



and 



previo usly claim ed. The broad correlation seen in the Gon-^ 
. Stanford et al. (2000| ) data is consistent 



or Seyfert) optical classifications, the remaining three are 
narrow line objects. Since the numbers of NVSS detected 
sources are low, it is difficult to determine accurate relations 
for individual source types. Nevertheless determinations of 
the q parameter, which in this paper is calculated from the 
rest-frame fluxes as 



logio{S(iQ^m./ Si 



(7) 



shows that the mean value for the NVSS detected 
HyLIGs is 1.66. This is somewhat lower than the typi- 
cal value determined for starburst galaxies and radio-quiet 
quasars of (~ 2.3). However, if the sample is split into the 
radio-loud TEXAS HyLIGs and the remainder we flnd that 
(a) the TEXAS HyLIGs have a median q value (0.89) which 
is consistent with the q ratios found for radio-loud quasars 
(0-1). (b) for the remaining HyLIGs with NVSS detections 
the mean q parameter is 1.94 which lies very slightly below 
the typical range for starburst and radio-quiet quasars of 
(2-2.6). However, this mean is much lower than the mean 
derived for LIGs an d UL IGs of q — 2.34 and q — 2.28 from 
Condon et al. (1991) and Stanford et al. (200Cl| ) respectively. 
This could be explained by the following postulates: 

• The NVSS-detected HyLIGs all contain AGN which 
contribute a significant fraction to the galaxy's total radio 
power, thus reducing the value of q. In ULIGs the radio 
power of any obscured AGN (if present) is minimal and 
hence radio emission from starburst regions is dominant. 

• HyLIGs have a skewed IMF with a higher proportion 
of high mass stars (which are prone to become supernovae) 
than ULIGs and therefore have a higher supernovae rate but 
maintain the same infrared power. Thus the radio emission 
(due to SNe rather than AGN) of HyLIGs is higher than 
that of ULIGs. Hence the value of q is lower in comparison. 

It is important to note that in Figure |^ using the ac- 
tual radio fluxes rather than NVSS upper limits will push 
the points to the left i.e. to higher values of q. Consid- 
ering the NVSS survey limit is 2.7mJy and the 60/im k- 
corrected fluxes together imply that the sources must have 
q values greater than 2 (with the only exception being IRAS 
F14218-t-3845). The calculated mean g- value is therefore bi- 
ased by HyLIGs with lower values of q. 

Therefore the postulates described above apply only 
to HyLIGs with lower q values than is expected for radio- 
quiet sources (i.e. less than 2.0). There are six such sources 
(excluding the TEXAS sources). It is possible that these 
sources are from a class of intermediate radio power AGN 
(i.e. in concordance with the first postulate). Despite the 
well known radio bi-modality of the quasar population, such 
sources have been shown to constitute a similar fraction of 
the quasar population as the relativistically boosted rad io- 
loud quasars ~ 10% ( |Falcke, Sherwood fc Patnaik 1996| )p^ 



and therefore is not unreasonable to expect such sources to 
be present within a HyLIG sample. 

Four HyLIGs have lower limits of q in excess of the 
range expec ted for ULIGs (i.e. > 2.6). This fraction is con- 



with that seen for HyLIGs. The majority of the HyLIGs 
which have NVSS detections have AGN (QSO, BALQSO 



sistent with Condon et al. (1991) who also found that the 
BGS LIGs and ULIGs had q values which were within the 
expected starburst range or higher. The remnants of Type 



not true for the new sub-millimetre HyLIGs. 
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Falcke et al. (199!:) attribute the radio-intermediate quasars 



as being due to relativistic boosting in RQQs. 
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Figure 6. The figur e above shows the correlation between the calculated 60^m and 1.41GHz power. The plot includes 40 ULIGs from 
Condon ot al. (1991). It also contains data for HyLIGs from the /50-HyLIG sample and from the compilation of known HyLIGs from 



RR2000. Plot (a) includes all sources from RR2000 and plots a radio upper limit of 2.7mjy if the source was not detected in the NVSS 
survey (2.7mjy is the completeness limit of the NVSS survey ). Plot (h) includes o nly those HyLIGs with confirmed NVS S counterparts. 
Overplotted are lines representing the median 9 values for the Condon et al. (1991) sample (qc), the Stanford et al. (2000) sample (qsr), 
the four radio-loud hyperluminous Texas sources (qth) and the remaining hyper luminous sources {qh)- 



II and Type lb SNe ac celerate the m ajority of relativistic 
electrons in starbursts (Condon 1992). Over the lifetime of 
the starburst, the number of relativistic electrons gradually 
increases with each new SNR [the lifetime of the relativistic 
electrons (10*j/r) is longer than t hat of the SNe (months- 



lO^yr) 



Condon (1992[ )]. Consid- 



years) and SNRs (~ 10^ 
ering this, it is plausible that the HyLIGs with the higher 
g- values are the younger starbursts which have had less time 
to build upon the relativistic electron reservoir. Therefore 
their radio emission is lower than that of older starbursts, 
causing the value of q to be higher. Further radio follow- 
up of the high-g HyLIGs will be performed to confirm their 
radio weakness. 



7 CONCLUSIONS 

Infrared Space Observatory measurements of hyperluminous 
infrared galaxies have been performed to ascertain the na- 
ture of these extreme infrared objects. 

a) Four objects, with faint IRAS detections, have been 
confirmed to be real infrared sources. 

b) The ISO observed sources constitute some of the 
most luminous sources in the Universe {Lir > 1 x 

c) The broad-band spectral energy distributions of these 
objects have been compared to starburst and Seyfert mod- 
els. The IR emission of IRAS F00235-hl024 and IRAS 
F14218-I-3845 are predominantly starburst fuelled whereas 
predominant AGN fuelling is seen for TXS0052-I-471 and 
IRAS F15307-^3252. 

d) The only radio loud HyLIG of this sample 
(TXS0052-I-471) has emission consistent with a dust torus 



model (EHY95) with an inclination angle of 26°. This ob- 
ject is of most interest since it is part of a subset of faint 
IRAS F sources which were correlated with the TEXAS Ra- 
dio Survey to identify possible extreme sources. Prior to the 
detections by ISO, the reality of TXS0052-)-471 had not been 
confirmed. 

e) Photometry data of sufficient quality has been ob- 
tained for all four HyLIGs to develop more accurate emission 
models for hyperluminous IR sources. 

f ) The radio- far infrared luminosity correlation has been 
verified to continue to previously un-investigated radio and 
infrared luminosity powers. The mean q value for the radio- 
quiet sources (1.94) is lower than that previously determined 
for ULIGs (2.34) and indicates higher radio luminosities for 
HyLIGs. 

g) All the HyLIGs requiring a starburst compo- 
nent to explain the IR SEDs have SFRs in excess of 
10'^'^ MQyr~^ h'^Q (adopting either a Salpeter or Miller-Scalo 
IMF). 

h) We postulate that most HyLIGs require a combi- 
nation of starburst and AGN components. The results of 
this study imply that better sampling of the IR emission of 
HyLIGs may reveal that an AGN component (not necessar- 
ily dominant) is required to explain all the emission from 
the NIR to the sub-millimetre. 
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Name 


z 


Spec. 


■560 


51.4 


Logxo(PiB) 


Logio(Pi.4) 


q 


logio{LiR/LQ) 






Type 


(mJy) 


(mJy) 










IRAS F00235+1024 


0.58 


NL 


428 


2.6 


27.084 


24.610 


2.475 


13.016" 


TXS0052+4710 


1.93 


QSO 


223.2 


150.5 


28.223 


27.445 


0.778 


13.834" 


SMMJ02399-0136 


2.803 


S2 


428 


< 2.7 


28.983 


26.029 


> 2.954 




IRAS P07380-2342 


0.292 


NL 


1170 


< 2.7 


26.798 


24.016 


> 2.781 




TXS0749+211 


2.2 


QSO 


155.4 


31.6 


28.231 


26.883 


1.348 


14.226 


IRAS F08279+5255 


3.91 


BALQSO 


511 


< 2.7 


29.496 


26.321 


> 3.175 




IRAS P09104+4109 


0.44 


S2 


525 


15.9 


26.875 


25.150 


1.725 


12.767 


iXbU9i5+Ui9 


2.57 


QSO 


158.7 


138.0 


O O A A r\ 

28.440 


27.661 


0.779 


14.545*' 


IRAS F10026+4949 


1.12 


SI 


266 


< 2.7 


27.631 


25.213 


> 2.418 


13.513*' 


TXSlOll+142 


1.55 


QSO 


225 


165.7 


27.953 


27.291 


0.661 


13.564" 


IRAS F10214+4724 


2.286 


S2 


190 


< 2.7 


28.368 


25.849 


> 2.519 


14.137" 


PG1148+549 


0.969 


QSO 


196 


4.3 


27.327 


25.286 


2.041 




PG1206+459 


1.158 


QSO 


463 


< 2.7 


27.912 


25.243 


> 2.668 




IRAS F12207+0939 


0.68 


QSO 


180 


5.8 


26.885 


25.100 


1.785 


13.711 


IRAS F12358+1807 


0.26 


BALQSO 


262 


5.6 


26.031 


24.230 


1.801 


12.057*' 


PG1248+401 


1.03 


QSO 


224 


< 2.7 


27.457 


25.139 


> 2.319 




IRAS F12509+3122 


0.78 


QSO 


218 


< 2.7 


27.123 


24.890 


> 2.233 


13.144 


IRAS F12514+1027 


0.3 


S2 


712 


8.5 


26.609 


24.538 


2.071 




PG1254+047 


1.024 


QSO 


307 


< 2.7 


27.587 


25.133 


> 2.454 




IRAS F13279+3401 


0.36 


QSO 


1028 


< 2.7 


26.956 


24.202 


> 2.754 


12.920" 


IRAS P14026+4341 


0.324 


QSO 


609.8 


< 2.7 


26.621 


24.108 


> 2.512 




H1413+117 


2.546 


BALQSO 


230 


7.9 


28.589 


26.410 


2.179 




IRAS F14218+3845 


1.21 


QSO 


36.16* 


< 2.7 


26.857 


25.282 


> 1.575 


12.856" 


IRAS b 14481+4454 


0.66 


S2 


190 


12.1 


26.875 


25.393 


1.482 




IRAS F14537+1950 


0.64 


SB 


283 


< 2.7 


27.013 


24.714 


> 2.299 




IRAS F15307+3252 


0.93 


S2 


280 


8.1 


27.434 


25.524 


1.910 


13.255" 


FFJ1614+3234 


0.71 


NL 


174 


< 2.7 


26.918 


24.806 


> 2.112 




PG1634+706 


1.334 


QSO 


318 


2.4 


27.919 


25.318 


2.601 




IRAS P18216+6418 


0.3 


NL 


1128 


91.9 


26.809 


25.572 


1.237 




IRAS F23569-0341 


0.59 


NL 


347 


< 2.7 


27.012 


24.641 


> 2.371 


12.845" 



Table 3. Table of HyLIGs from RR2000 sample with redshift, spectroscopic type (taken from RR2000), the flux at 60/im SgOi the flux 
at 1.4GHz Si. 4, the power at GOfira logioPeofim, the power at 1.4GHz logioPi.iGH z and the q parameter where q = /o3io(S6o/>S'i.4). 
Finally, for the sources from the originally proposed /50-HyLIG sample, loo-in(Lrn/ L,^,) is given wher e " indicates an integration of 
a best-fitting SED over 42. 5-122. 5/im ' luminosity calculated as per the definition of 
60/^m fiux is extrapolated from the best-fitting SED model. 
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